Heat transfer between nanoparticles: Thermal conductance for near-field interactions 
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We analyze the heat transfer between two nanoparticles separated by a distance lying in the near- 
field domain in which energy interchange is due to Coulomb interactions. The thermal conductance 
is computed by assuming that the particles have charge distributions characterized by fluctuating 
multipole moments in equilibrium with heat baths at two different temperatures. This quantity fol- 
lows from the fluctuation-dissipation theorem (FDT) for the fluctuations of the multipolar moments. 
We compare the behavior of the conductance as a function of the distance between the particles 
with the result obtained by means of molecular dynamics simulations. The formalism proposed 
enables us to provide a comprehensive explanation of the marked growth of the conductance when 
decreasing the distance between the nanoparticles. 
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I. INTRODUCTION 

The study of energy transfer mechanisms at the 
nanoscale-"^? has aroused increasing interest due to the 
emergence of the interdisciplinary field of nanoscience 
where such wide-ranging fields as for example solid state 
physics^, nanothermodynamics^i^ or electrical engineer- 
ing^ coexist. One of the basic problems in this field is to 
determine the energy exchange between two nanoparti- 
cles (NPs) at different temperatures. The way in which 
this energy is transfered depends crucially on the dis- 
tance between the particles. For sufficiently large dis- 
tances, heat exchange proceeds via thermal radiation, 
through emission or absorption of photons whereas at 
smaller distances recent molecular dynamics simulations 
have shown that Coulomb interaction (near-field radia- 
tion) is the dominant mechanism^. 

For near-field interactions, the thermal conductance 
was calculated under the assumption that both NPs 
behave as effective dipoles at different temperatures^. 
Hence, since these dipoles undergo thermal fluctuations, 
the fluctuation-dissipation theorem (FDT)^'— liiii^ pro- 
vides the energy which dissipates into heat in each NP. ft 
was found that the heat and therefore the conductance 
varies according to a very different behavior from 
the one observed in the case of thermal radiation: d^^ . 
Molecular dynamics simulations agree with the dipole- 
dipole model when the two NPs are separated by a dis- 
tance on the order of a few nanometers. However, near 
contact the conductance deviates dramatically from the 
prediction of the dipole model, as the simulations show. 
This behavior is a consequence of the fact that when 
particles become very close the position of the atoms are 
highly correlated, consequently the charge distributions 
become nonsymmetric and cannot be described merely 
as two interacting dipoles. To account for this distor- 
tion of the distribution of charges a more general formal- 
ism which focuses more convoluted interactions involving 



higher order multipoles aside from the dipoles is required. 

Our purpose in this paper is to provide this general 
formalism enabling us to analyze the behavior of the con- 
ductance beyond the dipolar approximation. We will use 
the linear response theory to derive an expression of the 
FDT for the fluctuations of the higher order multipoles. 
In particular, we will focus on the quadrupolar contribu- 
tions to the conductance which are able to reproduce the 
behavior observed in the simulations for some sizes of the 
NPs. 

The paper is organized as follows. In Section 2, we 
present the multipolar expansion of the Coulomb forces^^ 
between both NPs and derive a general expression of the 
FDT valid for multipoles of any order which leads to the 
heat transfer between the NPs. In Section 3, we analyze 
the particular case of quadrupolar contributions and de- 
rive the expresion of the conductance. We compare our 
result with the molecular dynamics simulations^. Finally, 
in Section 4, we emphasize our main conclusions. 



II. HEAT TRANSFER BETWEEN TWO 
NANOPARTICLES 

In this section, we will study the near-filed radia- 
tive heat transfer flux between two NPs which interact 
through Coulomb forces. 



A. Multipolar expansion 

To analyze the Coulomb interaction between two NPs 
(see Fig.[T]) it is necessary to know the charge distribution 
inside each of them. This can be performed by specifying 
their multipole moments so that the multipole moment 

can be defined asi^ 



of order n of the NPi, m|")^ 
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FIG. 1: Schematic illustration of the interaction between two 
NPs (NPl and NP2) with temperatures Ti and T2, respec- 
tively. Each NP is assimilated to a multipole moment (mo- 
ments Mj'' and Mj™'), and are separated by a distance d 
between their centers. 



where Cr is the charge at the position r inside the NP 
and the ^^"■'(r) are symmetric irreducible tensors (see 
the Appendix I for more details) where a = (ai, . . . , Q!„) 
and Uj = 1,2,3 for j — Thus, the case n — Q 

corresponds to the monopole mI? (r) = X^r'^r, n ^ 1 



is related to the dipole moment M'l^y^{r) = X^r^rJ^ai, 
and n = 2 for the quadrupole moment M^V^.^(r) = 
1/2 Br (3rQjQ2 — 5aicx2)- Hcncc, in terms of the spher- 
ical surface tensors Ya^\v), given through Eq. (jA-3[) . 
and the unit vector f related to r, Eq. ^ adopts the 
equivalent form 

^(SoW = ;j!E^rr"yi")(f). (2) 

r 

The above mentioned interaction between these NPs 
modifies their respective Hamiltonians. The interaction 
between NPz with NPj introduces a time-dependent per- 
turbation H(AA\ in its Hamiltonian which can be written 



Hi]) 

as a multipolar expansioni^: 



m=0 



(m) 



(3) 



with c„i = 1/(2to — 1)!! and stands for the full con- 

W (1,3) W;a (jj);q 

y{m) 



traction of indexes, M^") V^^-l = m/T^ y/",^ . In 
addition 



(4) 



with V(j (d) being the interaction potential between 
both NPs and d the separation between their centers. 
In terms of the first contributions, the perturbation can 
be expressed as 



1. 



(2) ^ <ri^) 



(5) 



where — v|^^^ is the electric field induced in the NPi, 

— V^^ is the gradient of this induced field, and M^^^ is 
the conjugated quadrupolar moment. 



Likewise, the electrostatic potential admits a multipo- 
lar expansion as well 



(6) 



which expresses the fact that the potential acting on the 
NPi depends on the charge distribution in the NPj. Here, 
m|Jj^ are the multipolar moments of the NPj and 



1 



— - — r^")(d) 



(7) 



is the Green propagator, with d the vector connection 
the centers of the particles and d the corresponding unit 
vector. Thus, from Eq. ^ 



n=0 



where G^™a"' is defined through 



(8) 



(9) 



B. Heat transfer from the fluctuation-dissipation 
theorem 



In the linear response regime, the multipolar moments 
can be expressed as 

oc 

where p("'™^ (w) are the multipolar polarizabilities which 
may in general depend on frequency. 

The energy transferred between the particles and con- 
verted into heat can be obtained from the linear response 
theorjiiiii^. One obtains (see the Appendix II) 



n.m— 



(11) 



where the symbol * stands for the complex conjugated, 
and the brakets express thermal average. 

According to Eq. ^ , the term in Eq. (|lip containing 
the thermal average can be transformed as 



00 00 



n,m— 



m[P)* gC'") pj;;"^ G("+*'-) Mf ) 
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l,k=0 



(12) 
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where we have defined 

oo 

= E Gi^''"^ ^\T^ ® . (13) 

n,m— 



Moreover, from Eqs. (|A-3[) . ([7]) and ([5]) one can prove 
that the s|^.'^'^'' are symmetric tensors. Therefore, making 
use of Eq. Eq. HIl) becomes 



-ilusq 



E 



ClCk X 



{(M«;0S|;f 0MW)-c.c.}. (14) 



The dependence of the energy transferred on the distance 
d resides in \ as follows from Eq. (fT5|) and the ex- 
pression of the propagators given through Eqs. ([7]) and 
^ . The multipole-multipole correlation can be obtained 
by using the FDT^^^^^ 



TTUJClCk 



(15) 

where e(w,T,) = ftw {1/2 + 1/ exp(fiij/fcTj - 1)} is the 
mean energy of an oscillator. As an illustration, for the 
dipolar case^, we obtain 



,(1,1) 



{^'^'-^'if*)Q{^.T,\ (16) 



where P^^^' is the dipole-dipole polarizability which we 
assume to be given through 

45£^ = «wHaS (17) 

and A^^^j, = (^q^^j, with a^^i^) = a(i)(w) + «a'(i)(^)- 
Hence, Eq. ([TB)) becomes 



«!>4)l/5> = ^«(.) To^.,,,. (18) 

For the quadrupolar case one has 



M(^)*M(^n = ^ fpL^f ^ - P^''* ) e(c.,T.), (19) 



ePg^^i. 
given through 



(2 2"). 

where P^^-j' is the quadrupole-quadrupole polarizability 



with /3(,)(w) = /3[j)(t^) + j/3('i)(w) and 

(2) 1 



(20) 



(21) 



is the isotropic tetradric. Thus, with Eqs. ([20)1 and (|2ip . 
Eq. (fT9| is written as 



(i);a (i);/3 



2eo 



,/3[',)(c.)e(c.,T,)Ai^)^. (22) 



7rw(c2)- 



Up to the quadrupolar order one has to take into account 
also the cross correlation dipole-quadrupole 



(i)*A/r(2)\ _ -'gp 



(0 (») 



fp(i,: 



,(1,2)^ 

(0 



e(c^,T,; 



(23) 

fl 2") ■ 

where P^^^' is the dipole-quadrupole polarizability, given 
through 



«';2;3 =7w(^)na/,/3i,/3: 



(1) 



with 7(,)(w) = 7|,)(^) + *7(,)(^) and 



□ 



(1) 



(24) 



(25) 



an isotropic skew-symmetric tensor. From Eqs.(P ^ - (P5)) 
it follows 



(i);a ('i):/3 



2e, 



7[',)(c.)e(L.,T,)e«„0„&. (26) 



It must be emphasize that the FDT, Eq. p^ . applies 
whenever the charge distribution of each particle in the 
presence of mutual interactions has reached equilibrium 
with the heat bath. When this is not the case, in the non- 
aged regime^^, one can still use a similar expression of the 
FDT in terms of an effective temperature. This can be 
done through a generalized Langevin equationi^, which 
takes into account the heat exchange between the NPi 
and its thermal bath. Relating the momentum variance 
of the NPi with its temperature by the equipartition the- 
orem, one can obtain the effective temperature (2^i^/(i)) 
through the response of the system due to fluctuations 
of the multipolar momentsi^. Thus, 



T 



{i,k} 

ffii) 



To + {To - n] 



, (27) 



where Tq stands for the initial temperature of the NPi 
and Tj, is the bath temperature. This expression shows 
that when the multipole moments of the particles are un- 
correlated, i.e. when both particles equilibrate indepen- 
dently at two different temperatures, the effective tem- 
perature coincides with that of the bath. This is the 
situation addressed in this paper. 

The effective temperature, defined as that for which 
the system would equilibrate, is a parameter measuring 
the distance to the stationary state in which both par- 
ticles reach two different temperatures. It can be also 
calculated using a relaxation mode l^^'^^ . 
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III. THE THERMAL CONDUCTANCE 

In this section, we will calculate the thermal conduc- 
tance between the two NPs in the presence of quadrupo- 
lar contributions. To this end, we start by writing Eq. 
(fT4| as follows 

— ILO£0 



Lk 



Qi—*j 



l,k=0 {13} ,{a} 



(28) 



where from Eqs. dH) and (fT3 



,(Z,fc) 

'a);/3,. 



(47r£o) 



^ OO 



7i,m— 



with 

{7} M 

yj^'^Hd). (30) 

Hence, by substituting Eq. (|5D|) into Eq. ([^ and the 
resulting equation into Eq. (|28p we obtain 

00 00 

^ ' Lk=0 n.m=0 



(29) 



(31) 



where 



{Lk) 



(32) 



Therefore, from Eq. PT|) we obtain the net heat flux 
between both NPs 

Q12 — Qi^2 — Q2^1 

00 00 

= E -'-^ E X 

^ ' Lk=0 n,m=a 



.2) (2->l) 



(33) 



In view of Eq. (j25p and the symmetric character of 
the spherical surface tensors given through Eq. (jA-3P 
one can prove from Eq. pO|) that 



= (34) 

Moreover, it can be shown that when m + n = 2p + 1 
{p > n), pIJ)'™^ is proportional to an isotropic skew- 
symmetric tensor D^p^ of order 2p+ 1 which satisfiesi^ 

□^"1 , = -D^P^ , (35) 



Therefore, by symmetry reasons only coefficients C^^^^-^ 
for which n + m — 2q and I + k = 2s, with g and s two 
positive integers, contribute to the heat flux. Hence, up 
to quadrupolar order we can write from Eq. (|3ip 



1 



_ f^l, 1.1, 1^-6 , /^l, 1.2,2 

+9Clf_:]^l)d-' + Cl^fd-'° + ...}. (36) 



Therefore, from Eqs. ([36]), and (|A^ - (fA^ we 

arrive at 



3 

in'' 
15 



140/3^'^) (a;)/3;'^.)(u;)ci-i° + . . .} e(^,r,).(37) 
and consequently 



;i2 



(uj) = Ql^2 - Q2^1 



8^ {«(1)"(2)^ ' + 



Ae 



45 
8^ 



{("(2)/3(i) +3a(i)/3(2)) 0(w,ri)- 



(a'(',)/3^2) + 3a'('2)/3^;)) e(c^, T2)} rf-^ (38) 



where AG = {e(w,Ti) - e(tJ,T2)}. 

When NPs are at the same temperature T, Eq. 
reduces to 



^121 



,(^) = — (a'(',)/3^2) - a'l,/l,^) d-^e{uj,T), (39) 
whence since the system is in thermal equihbrium 

QuHduj = 0. (40) 



In the general case, i.e. out of equilibrium, we can lin- 
earize Eq. ([55)1 with respect to the temperature diference 
AT = Ti ~ T2 in order to obtain the conductance given 
through Gi2(To) = dQu/d AT\j,^^j.^^j.^. We obtain 

Gi2(ro) = e'{L,,To)[al^a'l^^d-' + 

60 (a(i)/3(2)+«(2)/5(i))'^"' + 
5/3[;)(c.)/3^2)(c.)d-i"}dc., (41) 

where Tq = (Ti-f T2)/2 is the average temperature, which 
corresponds to the final equilibrium temperature that 
two bodies would reach when brought into contact and a 
heat flow established between them^^. 

In the expression we have obtained for the conduc- 
tance, we can identify the following contributions: 
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FIG. 2: Thermal conductance G12 vs distance d reproducing 
the molecular dynamics data obtained by Domingues et al^. 
The grey points represent the conductance when the particles 
with effective radius 0.72, 1.10, and 1.79 nanometers are in 
contact. The dotted lines show the analytical result obtained 
by these authors. The values of the conductance at short 
distances have been given through the grey dash-dotted lines. 
The inset shows different guide behaviors as a function of the 
distance having different values of the exponents. 



(i) Dipolar 

Gil' (To) = ^ (^°° e'(c^, n)al,^^al,^<h?j d^'. (42) 

which coincides with the expression obtained in Ref<^. 

(ii) Quadrupolar 

1 f°° r / 

«(2)/^(i)) + 

^Pl^^iu;)P'l^^ic.)d-'"}cL;. (43) 



In order to verify our results, in Fig. [5] we reex- 
hibit a graph obtained by Domingues et al^ extending 
the logarithmic scale for conductance in the more usual 
form. This graph displays the thermal conductance as a 
function of distance between the NPs, both with radius 
i?, in three significant situations: in mechanical contact 
(d = 2i?), in the intermediate region shortly before con- 
tact (2i? < d < 4i?), and in the most distant region 
{d ^ 4i?) where the near-field interaction is still valid. In 
this situation, the results corresponding to the grey dot- 
ted lines show the behavior d~^ which was obtained in 
RefjS'. Our results are in broad agreement for this region 
where the dipolar domain is present. When the parti- 
cles are close togheter, their charge distributions becomes 



very disorderly and higher orders than dipolar interac- 
tions come into play in the calculation of the thermal con- 
ductance. In this case, as predicted by Domingues et al4 
the thermal conductance is about 4 orders of magnitude 
larger than that of the dipole model given in Eq. (|42p . 
In more extreme conditions when the particles come into 
contact to each other, the same authors also predicted 
that the conductance would be 2 to 3 orders of mag- 
nitude lower than the conductance just before contact. 
These numerical predictions are covered by the result we 
give in Eq. (|4ip where one can see that the dominant 
contribution [d^^^) is 4 orders of magnitude lower than 
the dipolar case [d^^) while an intermediate case would 
give a value d~^. 

It must be stressed that we have obtained the conduc- 
tance up to quadrupolar order, nonetheless through our 
formalism it is possible to obtain the conductance for any 
order of multipolar interaction. 



IV. CONCLUSIONS 

In this paper, we have presented a theory to explain 
the exchange of energy between two NPs at different tem- 
perature. Our theory provides a general formalism based 
on the multipolar expansion of the electrostatic field in 
order to study heat transfer between two NPs for arbi- 
trary small distances provided that the FDT be satisfied. 
However, out of the FDT regime and when the system 
possesses fast and slow degrees of freedom it is possible 
to formulate a FDT in terms of a non stationary effec- 
tive temperature which depends on the slow degrees of 
freedor 



-16.17 



We have found that our analysis of the heat inter- 
changed between two NPs separated by a few submicrons 
agrees with the explains the rapid growth of the conduc- 
tance observed in the simulation^, even when the NPs 
are in contact. Hence, we are able to provide a compre- 
hensive explanation of the numerical results reported in 
Ref^. 

The formalism presented could also be applied to other 
situations such as the radiative heat transfer between a 
small dielectric particle and a surface^^ and the study of 
the optical forces due the radiation of a thermal source^S, 
enabling us to go beyond the dipolar approximation. 
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APPENDIX I 



The remaining coefBcients are obtained in similar way 



In this Appendix we present some of the properties of 
the spherical surface tensors yi"-* (f ) related to the n-rank 
Cartesian tensor 



X(")(r)-(-ir 



gn 2 



daa^ ■ ■ ■ daa^ r 



(A-1) 



introduced in Eq. ([T]). The tensor x!^\r) are given in 
terms of the unit vector f related to r as^ - 

Xi")(r) =r-("+i)yi")(f). (A-2) 

This spherical surface tensors can be expressed as 

ri")(f) = (2n - 1)!! ...f„„, (A-3) 



where "Tai ■ ■ ■ ^re the symmetric irreducible tensor 
constructed with the components of f . The first three 
symmetric irreducible tensors are 



(A-4) 
(A-5) 



(A-6) 



The spherical surface tensors satisfy the following prop- 
ertyi^ 

^yi")(f)yi")(f) = (2n-l)!!n!. (A-7) 

{a} 

By using Eq. ((30)) we can obtain the coefficients de- 
fined in Eq. From Eqs. (|T1])-(I121), §3), and 
one has 



C, 



1,1,1,1 _ -li^TT f ^(1,1) 



£0 



{13} {"} 



K>S.)-c.c.} 



m M 



(A- 



and 



^2,2,2 2 ^ ^ I ^(2 2) y y (2,2, 



{/3} {a} 



4/3(,)(-)/3(,) (-) E E ^0 

{/3} M 

4(7!!)(4!)/3;'^)(c.)/3^;.)(c.)e(c.,r,). (A-9) 



1,1,2,2 _ -^^TT I (14 



{/3} 



4a;„(.)/);.,(„)i:i:y<-"i^;f' 

{/3} 

360a'^^.)(a;)/3[V)(cc>)e(cc>,r,) (A-10) 



and 



C, 



= 360a'^^)(^)/3^.)(c.)e(c.,T,-). (A-11) 



APPENDIX II 

This Appendix is devoted to the derivation of the ex- 
pression of the energy dissipated corresponding to Eq. 
(fTTj) . In the adiabatic case, for a perturbation of the 
form 



(B-1) 



where Xj is a generalized displacement and fj (t) is a gen- 
eralized force, the change in the energy of the system is 
equal to the mean value of the partial derivative of the 
Hamiltonian with respect to time. Since only the pertur- 
bation H in the Hamiltonian depends explicitely on time 
and dynamical observable of the system which is 

independent of time, we have 



dE/dt — —CjXjdfj/dt. 



(B-2) 



In the framework of linear response theory one assumes 
that 

/>oo 

CjXj{t)^ / ajk{T)fk{t - T)dT, (B-3) 
Jo 

a relation similar to Eq. pUj) . After introducing the 
Fourier transforms and combining Eqs. (|B-3|1 and (|B-2[) 
we can write 



1 



dE 

It " (27r)2 

{iu;')fj{uj')ajk{i^)fk{oj) 



dui / duj' exp{—i{uj + uj')t) x 



(B-4) 



Here, if the perturbation f acts over a finite time, the 
total energy dissipated is 



dE 2 

dt-j^ = — I du! I duj'S{uj + uj'){iuj') x 



fj{uj')ajk{cj)fk{uj) 



1 

2^ 



duj{-~iuj) f j{—uj)ajk{u}) X 



(B-5) 
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Since the total heat must be a real quantity which after performing the thermal average leads to the 

equation equivalent to Eq. (jlip 

duj{~iuj)fj{~uj)ajk{uj)fk{uj) = 



2tt 



J duj{f*ajkfk~ fja*kfk)ui. (B-6) iu 



Therefore, the heat at the frequency lo is given through 



Qiu;) = -- {{f;a,kfk) - ifkal^f*)) ■ (B- 



Qi^) = -Z ifi^^kfk - fkal^f*) , (B-7) 
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